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I . INTRODUCTION 


Some  missiles  have  warheads  which  carry  a large  number  of  bomblets 
as  payload.  These  are  deployed  at  a predetermined  altitude  and  disperse 
over  a wide  area.  The  bomblets  are  designed  to  have  self-induced  spin 
which  generates  a Magnus-type  lift;  they  are  often  called  Magnus  rotors. 
The  spin  also  arms  the  fuze  so  that  the  bomblet  is  prepared  to  explode 
on  impact.  A variety  of  factors  could  interfere  with  the  arming  process. 
Among  these  is  the  process  by  which  the  bomblets  are  separated  from  the 
warhead  and  propelled  into  the  free  stream  where  they  acquire  their  full 
spin.  This  process,  called  bomblet  ejection,  establishes  the  initial 
conditions  for  the  subsequent  flight. 

The  specific  missile  and  warhead  considered  in  analyzing  ejection 
is  the  Lance  with  its  M251  warhead.  It  carries  about  800  bomblets  which 
are  essentially  spheres  (typically,  diameter  = 0.0572m)  with  ridges  to 
impart  spin.  They  are  deployed  at  low  supersonic  or  high  subsonic 
missile  speeds  after  line  charges  cut  the  missile  skin  into  two  panels. 

In  older  designs  an  active  ejection  system  was  used;  e.g.,  "slings," 
activated  by  the  line  charges,  propelled  the  bomblets  away  from  the  mis- 
sile. In  more  recent  designs  an  active  ejector  is  not  employed  but  the 
bomblets  are  usually  ejected  successfully.  However,  the  governing 
physical  process  was  not  understood. 

The  purpose  of  this  paper  is  to  construct  a model  of  the  bomblet 
ejection  process.  A model  of  the  process  is  constructed,  rather  than  a 
physical  description  from  first  principles,  because  there  are  not  enough 
data  for  the  latter.  It  is  specifically  intended  for  small  times;  i.e., 
when  the  radial  displacement  of  the  bomblets  is  of  the  order  of  10  d, 
where  d is  the  bomblet  diameter,  or  1 missile  diameter.  The  model  has 
proven  useful  in  the  explanation  of  certain  anomalies  in  the  performance 
of  the  warhead. 

Several  possible  mechanisms  for  ejection  are  examined:  inertia 
forces,  forces  due  to  explosive  line  charges,  pitching  motion  of  the 
missile,  and  the  roll  rate  of  the  missile.  The  last  of  these  is  found 
to  be  dominant.  The  important  elements  of  the  model  are:  the  flow 
field  after  the  panels  have  cleared  the  missile,  the  restraining  force 
on  the  bomblet  (which  introduces  the  dynamic  pressure)  as  it  crosses  a 
shear  layer,  and  the  roll  rate  of  the  missile.  A straightforward  calcu- 
lation gives  the  ejection  time  as  a function  of  roll  rate,  flow  field, 
restraining  force,  and  bomblet  position.  Despite  a number  of  simplifi- 
cations, it  is  felt  that  the  essential  features  of  the  ejection  process 
are  modeled  correctly. 

An  important  conclusion  is  that  the  time  for  ejection  is  a sensi- 
tive function  of  roll  rate  for  the  range  of  rates  used  in  practice.  The 
time  for  ejection  can  become  infinite  for  small  roll  rates,  which 
implies  that  the  bomblet  is  trapped  in  a pocket  of  separated  flow.  In 
reality  this  means  ejection  is  delayed;  the  bomblets  eventually  clear 
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the  missile,  at  times  so  large  that  the  model  no  longer  applies.  As  the 
roll  rate  approaches  zero,  the  model  predicts  large  ejection  times  for 
all  bomblets.  The  basis  for  the  model  no  longer  exists  since  the  effects 
which  were  negligible  must  become  dominant  if  the  roll  rate  goes  to 
zero.  An  alternate  model,  not  considered  here,  was  developed  which 
would  give  ejection  even  for  zero  roll  rate;  the  results  were  quite 
different  from  those  of  the  first  model. 

For  system  performance  studies,  a knowledge  of  bomblet  ejection  is 
useful  in  two  ways.  The  model  can  give  initial  conditions  for  a detailed 
study  of  bomblet  trajectories  and  spin.  Rather  than  a deterministic 
approach,  a statistical  one  can  be  used.  The  function  of  the  model  is 
then  to  provide  correlation  parameters;  this  will  be  discussed  in  a 
separate  report.  Its  success  provides  some  indirect,  partial  support 
for  the  validity  of  the  model. 

The  only  data  available  to  evaluate  the  model  are  the  time  of 
appearance  and  the  shape  of  the  distribution  of  the  bomblets,  as  deter- 
mined from  a few  film  records.  It  is  difficult  to  obtain  accurate 
readings  from  these  records  but  they  show  that,  typically,  the  bomblets 
first  appear  at  100  msec  at  which  time  they  are  about  0.61m  from  the 
missile.  Results  from  the  model  are  in  general  agreement  with  these 
observations.  In  one  case  the  film  records  were  clear  enough  to  measure 
the  shape  of  the  distribution  of  bomblets.  The  calculated  results  are 
consistent  with  the  measurement.  Thus  some  basis  for  confidence  in  the 
model  exists;  no  more  than  that  can  be  said,  at  present. 


II.  THE  BOMBLET  EJECTION  MODEL 

There  are  three  items  to  consider:  (1)  the  forces  acting  on  the 
bomblets  other  than  the  flow  field  forces,  (2)  the  flow  field,  and  (3) 
the  force  on  a bomblet  as  it  penetrates  the  shear  layer. 

The  packing  of  the  bomblets  in  the  warhead  must  be  discussed  first. 
From  the  manner  in  which  they  are  packed,  the  initial  distribution  of 
bomblets  will  vary  with  each  warhead.  Figure  1 shows  an  ideal  distri- 
bution which  gives  the  correct  total  number  of  bomblets;  3 rows  are 
shown  and  the  bomblets  are  numbered  for  later  reference.  In  practice, 
the  bomblets  are  not  necessarily  loaded  in  such  a regular  pattern  and 
distribution  would  be  somewhat  different,  but  unknown.  A choice  of 
initial  positions  must  be  made,  however,  and  the  one  shown  in  Figure  1 
is  convenient  for  calculations.  The  skin  and  nose  are  shown  by  dashed 
lines.  At  "event  time"  the  skin  panels  and  nose  are  blown  off,  leaving 
the  configuration  shown  by  the  solid  lines. 

The  missile  flies  at  small  angle  of  attack  with  a nominal  roll  rate 
of  5 rev/sec.  Only  supersonic  flight  at  event  time  is  considered.  It 
takes  about  10  msec  for  the  panels  to  lift  off  and  about  50  msec  to 
complete  the  separation  process.  The  fuze  ring  and  pedestal  may  or  may 
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not  be  lost  at  early  times.  After  the  panels  and  nose  are  removed  a new 
flow  field  is  established.  From  the  correlations  given  in  Reference  1, 
about  10  msec  are  required  to  establish  the  flow  with  the  fuze  ring  in 
place  and  probably  about  the  same  without  the  fuze  ring.  A typical  time 
of  10  msec  will  be  used  in  the  discussion  of  the  forces  that  were  exam- 
ined vis  a vis  ejection.  The  parameters  for  a specific  Lance  flight 
will  be  used  in  making  the  following  estimates. 

A.  Forces  on  the  Bomblets,  Excluding  Flow  Field  Forces 

Immediately  after  skin  release  the  forces  are: 

(1)  An  inertia  force  which  tends  to  move  the  bomblets  to  the  front. 
This  is  caused  by  the  sudden  change  in  drag  coefficient,  AC^,  because 

the  low  drag  nose  is  ejected.  An  adequate  approximation  is  AC^  = 1.0. 

The  axial  acceleration  of  the  bomblets  with  respect  to  the  missile  is 
then  47.5m/sec2.  In  10  msec  the  bomblets  are  displaced  axially  0.0023m. 
This  effect  will  be  neglected, 

(2)  Forces  due  to  the  line  charge.  Blast  and  detonation  waves 
cause  an  impulsive  force.  The  pressures  are  several  MPa!s  but  act  only 
near  the  charge  and  last  for  about  10  ysec.  This  can  be  neglected. 

From  the  burning  of  detonation  products  there  will  be  a residual  pres- 
sure in  the  warhead  volume,  estimated  to  be  2.07MPa.  This  is  relieved 
by  an  expansion  wave  propagating  in  from  the  moving  panels;  while  the 
wave  traverses  a bomblet,  an  outward  radial  force  acts  on  it.  A rough 
model  of  this  effect  gives  a maximum  acceleration  of  97.5m/secz  in  a 
radial  direction.  In  10  msec  the  bomblet  would  move  0.0048m,  if  the 
acceleration  were  constant.  But  it  lasts  for  about  1 msec.  Therefore, 
this  will  be  neglected. 

(3)  Pitching  motion  of  the  missile.  When  the  nose  and  skin  panels 
are  removed,  forces  act  on  the  missile  that  cause  pitching  motion.  With 
respect  to  the  missile  some  bomblets  will  accelerate  in  a direction 
normal  to  the  missile  axis.  The  information  necessary  to  calculate  the 
magnitude  of  this  effect  was  not  available.  However,  some  estimates  of 
it  were  made.  In  the  first  quarter-cycle,  which  lasts  about  170  msec, 
bomblets  on  the  windward  side  are  displaced  radially  whereas  those  on 
the  leeward  side  are  not  until  the  second  quarter-cycle.  Photographs  of 
a number  of  flights,  including  subsonic  and  supersonic  cases,  show  that 
the  time  for  the  bomblets  to  move  about  1 missile  diameter  (about  0.61m 
for  the  Lance)  ranges  from  80  to  130  msec.  Based  on  these  times  this 
mechanism  does  not  seem  plausible;  nor  does  it  from  an  estimate  of 
displacements.  If  no  other  forces  act,  in  10  msec  bomblet  la,  Figure  1, 
has  a significant  displacement,  about  0.02m,  but  bomblets  28a  or  7c 


1.  M.  S.  Holden 3 " Establishment  Time  of  Laminar  Separated  Flows 31  11  AIM 

Journal 3 Vol.  95  No.  11 3 November  1971 > pp.  2296-2298 . 
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would  be  displaced  0.002m,  which  is  negligible.  This  implies  that  the 
aggregate  of  bomblets  tips  during  ejection;  there  is  no  evidence  for 
this.  Although  it  has  not  been  established  that  this  mechanism  can  be 
neglected,  it  is  much  less  plausible  than  the  roll  rate  effect  and, 
therefore,  will  not  be  considered  further. 

(4)  Release  of  centrifugal  constraint.  The  nominal  roll  rate  is  5 
rev/sec.  When  the  skin  is  removed  the  bomblets  will  move  tangentially, 
in  a straight  line  at  constant  velocity  if  there  are  no  other  forces 
acting  on  the  bomblet;  see  Figure  2a.  For  the  inner  row,  a in  Figure  1, 
with  the  diameter  of  the  central  tube  0.178m,  the  initial  acceleration 
is  116m/sec2  (11.8  g) , the  initial  velocity  is  3.69m/sec  and  in  10 
msec  a bomblet  would  be  displaced  radially  0.0056m.  The  same  quanti- 
ties for  the  outer  row,  c,  are:  229m/sec2  (23.3  g) , 7.28m/sec,  and 
0.011m.  For  the  case  of  no  forces  the  radial  motion  is  obtained  from 
a simple  calculation.  The  distance,  r(t),  from  the  missile  axis  to  the 
bomblet  center,  at  time  t,  is  given  by 

r = a/cos  (tan  ^ 2irS3t)  (1) 

where  a = r(0)  and  Q,  is  the  roll  rate  in  rev/ sec.  Actually,  there  are 
forces,  described  in  Section  B below,  that  may  even  prevent  the  bomblet 
from  being  ejected  at  early  times.  These  forces  are  an  important  part 
of  the  model,  but  if  they  are  neglected  for  the  moment,  (1)  gives,  for 
= 5 rev/sec,  radial  positions  at  t = 100  msec  of  0.396m,  0.579m,  and 
0.762m  for  rows  a,  b and  c respectively.  These  radial  distances  are  in 
the  neighborhood  of  those  observed.  Of  course  the  bomblet  distribution 
will  not  be  correct  because  the  other  forces  have  been  neglected,  but 
even  this  overly  simplified  approach  gives  reasonable  results.  Release 
of  the  centrifugal  constraint  is  the  most  plausible  mechanism  for  bomb- 
let  ejection;  it  will  now  be  incorporated  into  a model. 

The  spin  of  the  bomblets  and  the  ridges  which  cause  it  are  ignored; 
they  are  treated  as  spheres  with  diameter  d = 0.0572m.  Only  the  motion 
of  the  center  of  mass  of  the  bomblets  will  be  described.  The  equations 
of  motion  are  written  in  a cylindrical  polar  coordinate  system  fixed 
with  respect  to  the  missile.  Since  the  missile  decelerates  this  is  not 
an  inertial  system.  However,  the  axial  deceleration  is  1/14  that  of  a 
bomblet  (or  W/C^A  is  14  times  greater)  and  will  be  neglected.  For  an 

inertial  system,  with  x the  axial  coordinate  along  the  center- line  of 
the  missile  and  r and  4>  polar  coordinates  in  the  plane  x = constant,  the 
equations  of  motion  are 


m (r  - r 4>2 ) = F^ 

(2) 

mr“  1 d(r2$)/dt  = F^  = 0 

(3) 

mx  = F 

(4) 
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with  initial  conditions,  at  t = 0,  when  the  flow  is  first  established. 


r = 

a 

o 

ii 

• u 

= 

0 

• 

<J>  = au> 

x = 

X 

o 

x = 0 

Here  m is  the  mass,  0.419kg,  of  the  bomblet,  and  a are  the  initial 

axial  and  radial  coordinates  of  the  bomblet,  co  is  the  roll  rate  in 
rad/sec,  and  the  applied  forces  are  F^,  F , F^.  Assuming  axial  sym- 
metry, there  is  no  force  in  the  4)  direction  and  (3)  can  be  integrated, 
giving 

r2ij>  = Cj  = a2u> 

This  equation  is  not  integrated  since  the  circumferential  motion  is  not 
needed  in  the  model.  It  is  used  to  simplify  (2).  Thus  the  equations  to 
be  integrated  are 


r = Ch  2/r3  + F /m 

1 r 

(5) 

x = F /m 

X 

(6) 

The  term  C^2/r3  is  the  apparent  centrifugal  force  which  accelerates  the 

bomblet  even  if  there  is  no  force;  i.e.,  if  F^  = 0.  In  that  case  (1) 

can  be  derived  from  (5);  but  this  is  the  hard  way  to  do  it.  If  F^  and 

F^  are  functions  of  r only  the  equations  can  be  integrated  by  quadrature 

after  using  the  substitutions  2r  = dr2/dr  and  x = rdx/dr.  In  general, 

F^  and  F^  are  functions  of  both  r and  x so  that  a numerical  integration 

is  required. 

B.  The  Flow  Field 

Although  some  parts  of  the  flow  field  will  certainly  be  unsteady  we 
shall  assume  steady  flow  here.  (There  are  some  analogies  with  flow  over 
a spiked-nose  projectile  which  are  useful  in  thinking  about  the  flow 
field.)  For  the  kind  of  bomblet  distribution  shown  in  Figure  1,  the 
steady  flow  near  the  bomblets  would  have  regions  of  attached  and  sepa- 
rated flow.  Some  schlieren  pictures  of  the  flow  over  a configuration 
similar  to  that  of  Figure  1 were  available.  Reference  2,  and  were  used 


2.  Private  corrmunioation  from  Mr .•  A.  Loeb. 
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as  a guide  for  approximating  the  flow  field.  The  shear  layer  indicated 

in  Figure  1 is  drawn  with  the  help  of  those  pictures.  The  shear  layer 

angle  with  respect  to  the  x-axis,  0 , varies  between  0 and  20°.  Actu- 

s 

ally,  slightly  negative  values  are  possible  and  the  pictures  of  Refer- 
ence 2 can  be  interpreted  to  give  0 =25°.  At  t = 0 most  of  the  bomb- 

lets  will  be  encased  in  the  shear  layer  but  the  flow  will  attach  on  some 
of  them;  these  will  be  subject  to  an  additional  force  not  considered  in 
the  next  paragraph*  Calculations  have  been  made  for  a range  of  values 
of  0^  and  for  a bomblet  on  which  the  flow  attaches.  The  shear  layer  is 

important  because,  until  they  penetrate  it,  the  bomblets  are  shielded  by 
it  from  the  flow  outside  the  shear  layer  which  has  high  dynamic  press- 
ure, q.  This  external  flow  could  be  obtained  from  a calculation  of  the 
inviscid  flow  over  the  fuze  ring  and  shear  layer.  Such  sophistication 
is  unwarranted;  the  flow  variables  immediately  outside  the  shear  layer 
were  obtained  from  the  flow  over  a flat-faced  cylinder.  Inside  the 
shear  layer,  q and  the  aerodynamic  forces  are  small.  In  Figure  2b  two 
typical  bomblets,  which  have  moved  from  their  original  position,  are 
shown.  The  only  significant  force  acting  on  them  is  the  centrifugal 
force  until  they  reach  the  shear  layer.  The  question  then  is:  can  a 
bomblet  penetrate  the  shear  layer? 

C . A Model  for  Penetration  of  the  Shear  Layer  by  a Bomblet 

The  complicated  interaction  between  a bomblet  and  the  shear  layer 
flow  must  be  modelled  in  a simple  way  that  retains  the  essence  of  the 
process.  Specifically,  an  approximation  to  the  forces  on  the  bomblet  is 
needed.  Three  assumptions  are  made  for  this  purpose.  First,  the  shear 
layer  is  replaced  by  a vortex  sheet;  i.e.,  an  inviscid  shear  layer. 
Second,  the  shear  layer  remains  fixed  while  a bomblet  penetrates  it. 
These  are  illustrated  in  Figure  2c;  the  small  velocity  inside  the  shear 
layer  is  taken  to  be  zero.  Let  the  components  of  the  force  on  the 
bomblet,  parallel  and  perpendicular  to  the  shear  layer,  be  F^  and  Fg; 

the  latter  is  called  the  restraining  force.  In  principle,  we  could 
calculate  the  force  on  the  bomblet  as  it  passes  through  the  vortex 
sheet,  treating  the  flow  as  inviscid  and  quasi-steady.  But  the  flow 
over  a segment  of  a bomblet,  position  2 in  Figure  2c,  is  three-dimen- 
sional and  such  an  elaborate  calculation  is  unjustified  and  unacceptable 
for  our  purpose.  There  is  one  position  of  the  bomblet  where  the  calcu- 
lation becomes  simple,  viz. , that  for  which  the  bomblet  center  is  on  the 
vortex  sheet,  position  3 in  Figure  2c.  Then  we  have  flow  over  a hemi- 
sphere. For  this  position  we  denote  F and  F by  F and  F . , respec- 

tively.  Experimentally  determined  pressure  distributions  over  a sphere 

can  be  used  to  calculate  F and  F u.  The  data  of  Reference  3 were  used 

sm  nh 


3.  F . D . Bennett , W.  C.  Carter  and  V . E.  Bergdolty  n Interferometric 
Analysis  of  Airflow  about  Projectiles  in  Free  Flight , " J . Applied 
Physics , Vol . 23,  No.  4,  April  1952y  pp.  453-469. 
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since  they  were  was  particularly  well  suited  for  this  purpose.  The  free 
stream  (or  missile)  Mach  number  is  Mm  = 1.82  for  a typical  flight  for 

which  results  will  be  presented.  Then  p = 0.5  pm  and  M = 2.13,  where  p 

and  M are  the  pressure  and  Mach  number  outside  the  vortex  sheet,  at  the 
position  of  bomblet  la.  Using  the  procedure  just  outlined  we  obtain 

F = 22. 2N 
sm 

F = 200N  . 

nh 

As  a check  on  the  calculation  note  that  2F  , = 400N  is  within  a few 

nh 

percent  of  the  drag  force  obtained  from  the  standard  for  a smooth 

sphere.  Because  of  the  ridges  on  the  bomblet  Fsm  and  F^  will  be  larger 

than  the  results  for  a smooth  sphere;  the  magnitude  of  the  increase 
depends  on  the  orientation  of  the  bomblet.  The  third  assumption  is  that 
the  forces  are  taken  to  be  linear  functions  of  s,  the  distance  of  the 
bomblet  across  the  shear  layer  and  that,  in  position  4,  Fs  = 0 and  F^  = 

drag  of  a sphere.  The  resulting  force  distributions  are  illustrated  in 
Figure  2d.  Some  obvious  refinements  to  these  force  distributions  could 
be  made  but,  because  of  uncertainties  in  the  physical  processes  we  are 
trying  to  analyze  and  in  other  elements  of  the  model,  this  is  not  justi- 
fied in  a first  approximation.  The  forces  that  enter  equations  (5)  and 

(6)  are  then 

Fr(s)  = -Fs(s)  cos  0s  + Fn(s)  sin  0s  (7) 

Fx(s)  = Fs(s)  sin  es  + Fn(s)  cos  es  (8) 

where 

s = (r  - r„)  cos  0 - (x  - x_)  sin  0 

v (r  s v 0 s 

and  r^  is  the  r coordinate  of  the  center  of  the  bomblet  when  it  is 

tangent  to  the  shear  layer,  as  in  position  1 in  Figure  2c.  Each  F in 

(7)  and  (8)  is  a linear  function  of  s,  or  constant;  they  are  zero  for  s 
< 0.  The  bomblet  accelerates  in  the  r direction  under  the  influence  of 
centrifugal  force  but  its  axial  coordinate  is  fixed  at  x = Xq  until  it 

becomes  tangent  to  the  shear  layer.  It  then  feels  the  effect  of  both  F^ 

and  F . If  its  radial  momentum  plus  the  centrifugal  effect  are  large 

enough  to  overcome  F , the  bomblet  penetrates  the  shear  layer  and  is 
ejected. 
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The  motion  is  easily  illustrated  for  the  special  case  6s  = 
because  (5)  can  be  integrated  once.  Consider  row  a so  that  r^  = a and 

s = r - a.  Then  from  Figure  2d,  where  s = d/2  at  position  3,  and  (7), 

t 

F = - F = - 2F  (r-a)/d  for  a < r < a + (d/2)  (9) 

r s sm 

= - 2F  [d  - (r-a)]/d  for  a + (d/2)  < r < a + d,  (10) 

sm  L 

Using  these  and  2r  = dr2/dr  in  (5)  yields 

r2  = a^w2  (a~2-r~^)  - 2F  (r-a)2/md  (11) 

v sm 

and 

r^  = a4a)2  (a“ 2-r~ 2)  + [d^  - 4d(r-a)  + (r-a)^]  F^/md  (12) 

for  the  two  ranges  of  r given  in  (9)  and  (10)  respectively.  For  the 
values  of  the ^parameters  a = 0.117m,  d = 0.0572m,  m = 0.419kg  and  using 

F = 35.6N,  r vs  r is  plotted  in  Figure  3 for  various  ft  = co/ 2tt ; t is  a 
sm  r 

parameter  along  each  curve  with  t = 0 at  r = a.  The  bomblet  penetrates 
the  shear  layer  for  ft  = 4.5  rev/sec,  curve  A,  but  does  not  for  any  of 

the  other  curves.  As  ft  decreases  the  bomblet  continues  to  escape,  but 

the  time  to  do  so  increases,  until  a limiting  value,  ft^,  is  reached; 

here  ft^  = 4.129  rev/sec.  For  this  curve  B,  (12)  has  a double  root  and 

the  time  to  reach  that  value  of  r is  infinite.  For  curves  C,  D,  and  E 

the  bomblet  would  oscilate  between  r = a and  the  other  for  which  r = 0, 

according  to  this  idealized  model. 

In  summary,  after  the  bomblets  have  been  exposed  and  the  flow 
established,  the  forces  are  defined  by  a knowledge  of  the  shear  layer, 
the  inviscid  flow  outside  the  shear  layer,  and  the  forces  acting  on  a 
hemisphere  at  the  local  flow  conditions.  The  magnitude  of  these  forces 
depends  on  the  dynamic  pressure.  Since  it  is  expected  that  this  model 
will  yield  only  estimates  or  be  used  for  comparative  purposes,  a precise 
knowledge  of  the  flow  field  is  not  necessary.  For  example,  if  the 
missile  is  at  angle  of  attack,  the  angle  and  position  of  the  shear  layer 
can  be  estimated  from  wind  tunnel  or  ballistic  range  shadowgraphs  of 
configurations  that  approximate  the  warhead  plus  bomblets.  In  any  case, 
to  obtain  useful  results,  the  parameters  will  have  to  be  varied  from, 
their  initial  estimates  because  of  the  number  of  simplifying  assumptions 
incorporated  in  the  model. 
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III.  RESULTS 


The  model  was  used  to  calculate  some  results  pertinent  to  ejection 
for  a typical  supersonic  flight  of  the  Lance.  The  velocity  at  event 
time  is  taken  as  609.6m/sec  (M^  = 1.82)  and  zero  angle  of  attack  is 

assumed.  The  input  quantities  are  xq,  a,  0^,  F , F ancl  r°H 
rate  in  rev/sec,  ft.  In  Figure  1 0^  varies  between  0°  and  20°.  To  sim- 
plify matters  the  results  are  presented  for  a constant  0^  for  all  bomb- 

lets  or  for  one  specific  bomblet.  It  is  more  tedious  to  do  the  calcula- 
tions using  the  0^  corresponding  to  each  bomblet.  Of  course,  as  the 

bomblets  penetrate  the  shear  layer,  0^  changes.  This  interaction  and 

that  of  one  bomblet  with  another  are  neglected;  the  motion  of  each 
bomblet  is  treated  without  regard  to  another.  Note  that,  for  a constant 
0^  = 5°,  the  shear  layer  would  be  a straight  line  joining  the  fuze  ring 

and  the  rear  wall  of  the  cargo  compartment. 

The  model  should  answer  the  questions:  will  the  bomblets  penetrate 
the  shear  layer  and  if  so  how  long  will  it  take?  The  ejection  time,  t^, 

is  defined  as  the  time  for  a bomblet  to  be  displaced  radially  one  diam- 
eter, across  the  shear  layer;  it  is  the  time  between  positions  1 and  4 

in  Figure  2c  and  is  also  indicated  in  Figure  4. 

For  bomblet  la,  see  Figure  1,  t vs.  ft,  for  various  0^,  is  presented 

in  Figure  4.  For  smaller  ft,  t^  is  a sensitive  function  of  0^.  In  fact 

each  curve  has  an  asymptote  for  some  ft;  i.e.,  t 00 , which  means  the 

bomblet  never  penetrates  the  shear  layer;  it  is  trapped  under  the  shear 
layer.  Figure  5 applies  to  all  bomblets  in  row  a and  shows  t^  vs  ft  for 

0 =0  for  a range  of  F above  and  below  the  value  of  22. 2N  derived  in 
s sm 

the  last  section  for  a smooth  sphere.  F = 0 implies  F^  = 0,  which  was 

discussed  in  the  last  section;  the  motion  is  determined  by  (1)  which 

also  gives  the  equation  of  the  dashed  curve  in  Figure  4: 


27rftt  = tan  [cos"1  ( a rj  ] 
e L a+d  J 


This  has  an  asymptote  at  ft  = 0.  For  each  F there  is  an  asymptote  at 

ft  = ft^;  the  relationship  between  ft^  and  F^  is  shown  in  Figure  6 for 

0^  = 0.  For  comparison,  the  values  for  two  other  0^  are  given:  for  0^ 

= 10°,  ft«  = 3.8  rev/sec  and  F = 23. 8N;  for  0 = 30°,  ft«  = 4.3  rev/sec 

t sm  s c 

and  Fsm  = 22. 2N.  The  asymptote  is  easily  explained  physically.  If  the 
radial  momentum  plus  the  centrifugal  effect  are  large  enough  to  overcome 
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the  radial  component  of  force,  the  bomblet  is  ejected;  i.e.,  t is 

0 

finite.  When  these  just  balance  tg  -*■  08 . Thus,  decreasing  ft,  increasing 

^sm’  0r  increasing  6S  will  force  the  bomblet  back  into  the  shear  layer 

and  it  will  not  be  ejected;  this  behavior  is  shown  in  computed  trajec- 
tories. The  correct  interpretation  of  this  result  is  that  the  ejection 
time  becomes  large,  by  which  time  the  angle  of  attack  of  the  missile  is 
so  large  that  the  model  can  no  longer  be  used  without  modification. 

Trapping  of  a bomblet  and  the  existence  of  ft^  was  explained  earlier 

using  Figure  5 for  the  special  case  0s  = 0.  In  principle  ft^  can  be 

found  by  determining  the  ft  for  which  (12)  has  double  roots.  Since  the 
roots  of  a quartic  must  be  found,  a convenient,  general  expression  for 

*-FsnP  can  ke  §i-ven;  however,  it  is  obvious  from  (12)  that  ft„  « p 

sm 

A convenient  expression  for  ft  (F  ) can  be  found,  where  ft  is  determined 
• sm  o 

from  r 0 at  r = a+d;  see  curve  C for  which  the  dashed  portion  has  no 
physical  meaning.  The  result  is 

^2  = (a+ci)  2 p //47f2a2m  (2a+d) 

o sm  v J 

Since  is  within  a few  percent  of  this  expression  is  adequate  con- 
sidering the  approximate  nature  of  the  model. 

The  sensitivity  of  t to  Q,  for  a given  F , is  not  the  same  for 

e sm 

all  bomblets,  even  for  a constant  9 . Figure  7 shows  t vs  ft  for  rows 

s e 

a,  b,  and  c of  Figure  1.  Because  0^  = 0,  t is  the  same  for  all  bomb- 

lets  in  a row.  The  results  for  row  a (already  given  in  Figure  5)  show 
that  t0  for  these  bomblets  is  more  sensitive  to  decreasing  Q than  it  is 

for  rows  b and  c.  This  is  because  the  centrifugal  force  acting  on  them, 
through  the  factor  C in  (5),  is  greater. 

The  definition  of  tg  does  not  account  for  the  time  it  takes  for  a 
bomblet  to  move  up  to  the  shear  layer,  t^,  which  can  be  obtained  from 

(1) . The  total  time  for  a bomblet  to  get  across  the  shear  layer  is 

ti  + t . This  is  shown  in  Figure  8 for  0 = 5°,  F = 22.2N  and  two 
x e s sm 

roll  rates,  3.8  and  5.0  rev/sec.  For  the  3.8  case  the  two  times  for 

bomblet  28a  are  t.  = 79.7  msec  and  t =26  msec;  for  bomblet  la.  t.  =0 

• e -1 
and  tg  = 66.6  msec.  t^  is  smaller  for  28a  because  of  its  greater  momen- 
tum when  it  reaches  the  shear  layer.  The  non-monotonic  behavior  for  row 
a with  ft  = 3.8  rev/sec  is  an  indication  that  the  first  few  bomblets  are 
close  to  the  asymptote. 
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A knowledge  of  t and  t^  are  useful  in  design  and  performance 

studies  and  the  existence  of  the  asymptote  is  an  important  result  from 
the  model.  Considering  the  various  appropriate  to  the  bomblets  we 

can  conclude  that  only  a fraction  of  the  bomblets  would  be  trapped  for 
certain  ft  and  F 

sm 


Trajectories  for  bomblets  la,  lb,  lc,  and  28a  are  shown  in  Figure  9 
with  ft  = 3.8  rev/sec  and  6g  = 5°;  the  velocities,  relative  to  the  mis- 
sile, at  t = 120  msec  are  also  listed.  Note  the  different  scales  in  the 
axial  and  radial  directions.  After  ejection  the  point  mass  trajectories 
are  computed  with  F^  = 0 and  F^  = 2F^  which  means  that  the  direction  of 

the  force  is  fixed  at  an  angle  6g  with  respect  to  the  missile  axis.  The 

lift  on  a bomblet  is  neglected.  The  three  dots  on  each  trajectory 
locate  t = 80,  100,  and  120  msec.  There  are  some  obvious  uses  for  the 
trajectories;  e.g.,  determining  if  a control  surface  is  struck  by  a 
bomblet. 

The  distribution  of  all  bomblets  at  t = 100  msec  is  shown  in  Figure 
10,  for  ft  = 3.8  rev/sec  and  Fgm  = 22. 2N.  The  shear  layer  angle  was  taken 

to  be  constant:  6 = 10°  was  chosen  to  give  the  best  agreement  with  a 

s 

test  result  discussed  below.  The  positions  of  rows  a,  b,  and  c are 
shown  by  curves;  the  first  and  last  bomblets  in  each  row  are  indicated 
by  a dot  and  are  numbered;  the  other  dots  show  the  positions  of  bomblets 
whose  number  is  a multiple  of  5.  The  spherical  bomblet  appears  to  be 
elliptical  because  of  the  different  x and  r scales.  The  distribution 
shows  that  some  bomblets  would  have  collided,  an  effect  which  is 
neglected.  For  a test  in  which  ft  = 3.8  rev/sec,  the  bomblets  could  be 
discerned  in  a photograph,  from  the  ground,  taken  at  t = 100  msec;  the 
bomblets  appear  as  a dark  cloud  near  the  missile.  The  leading  edge  of 
this  cloud  is  shown  in  Figure  10  as  a dashed  straight  line  connecting  the 
two  end  points  which  have  error  bars  attached.  The  quality  of  the 
photograph  did  not  permit  a more  accurate  measurement.  The  perspective 
of  the  photograph  does  not  allow  the  distribution  of  bomblets  to  be 
determined  but  the  leading  edge  can  be  compared  with  the  front  of  the 
calculated  distribution.  The  agreement  between  these  is  close  enough  to 
give  confidence  in  the  model. 

The  effect  on  the  distribution  of  changing  ft  is  shown  in  Figure  11. 
As  expected,  the  radial  displacement  of  the  bomblets  is  greater  for  the 
larger  ft;  the  axial  displacement  is  also  increased.  For  ft  < 3.8  rev/sec 
some  of  the  bomblets  would  not  appear  in  the  distribution  since  they 
would  have  been  trapped  by  the  shear  layer. 
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IV.  CONCLUSIONS 


The  problem  of  bomblet  ejection  from  missiles  which  do  not  have  an 
active  ejection  system  has  been  studied.  A model  was  formulated  for  the 
ejection  process,  the  crucial  elements  of  which  are  the  centrifugal 
effect  due  to  the  spin  of  the  missile  and  the  interaction  of  the  bomb- 
lets  with  the  shear  layer.  An  ejection  time  is  defined  which  can  be 
used  as  a measure  of  the  effect  of  various  parameters  on  the  ejection 
process.  For  a given  shear  layer  angle  and  restraining  force,  this  time 
becomes  infinite  as  the  roll  rate  decreases.  Taken  literally,  this 
implies  that  a bomblet  will  be  trapped  under  the  shear  layer.  Realistic- 
ally it  means  that  the  time  for  ejection  is  large,  at  which  time  the 
conditions  for  the  model's  validity  are  not  met.  Such  delayed  ejection 
can  affect  the  performance  of  the  system.  Bomblet  trajectories  and 
distributions  have  also  been  presented.  In  one  case,  the  calculated 
distribution  is  compared  with  the  shape  of  the  distribution  determined 
in  a test.  The  agreement  is  favorable  enough  to  lend  credence  to  the 
model . 
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Figure  1.  Sketch  of  a missile  warhead  containing  an  idealized  packing 
of  bomb  lets  in  three  rows:  a,  b,  and  c.  The  shear  layer 
is  formed  after  the  dashed  line  configuration  is  removed. 


V 


CONSTRAINT  REMOVED. 

BOMBLETS  MOVE  OFF 
TANGENTIALLY. 

V = tua  - 2 wfla 


Figure  2.  The  Elements  of  the  Bomblet  Ejection  Model. 

(a)  A cross-sectional  view  of  some  bomblets  on  the  central 
tube,  after  release  of  the  centrifugal  constraint. 


(b)  Two  typical  bomblets  moving  toward  the  shear  layer, 
with  only  centrifugal  force  acting. 
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Figure  2.  Continued 

(c)  Bomblet  penetration  of  the  vortex  sheet. 


Figure  2.  Concluded 
(d)  Assumed  force  distribution. 


21 


ro 


22 


Figure  3.  Phase  plane  diagram,  r2  vs  r,  of  the  motion  of  a bomblet  in  row 
for  es  = 0,  F = 35. 6N  (81bf),  m = 0.419kg  (.0287  slugs), 

a = 0.117m,  d = 0.0572m  and  n (rev/sec)  as  a parameter. 


110 

100 

90 

80 

70 

60 

50 

40 

30 

20 

Fic 


L I 1 

4 5 6 

£1  (rev/sec) 

4.  For  bomblet  la,  ejection  time  vs.  roll  rate,  with 
shear  layer  angle  as  the  parameter. 
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Figure  5.  For  all  bomblets  in  row  a,  ejection  time  vs.  roll 
rate,  with  restraining  force  as  parameter. 
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Figure  8.  The  total  time  for  each  bomblet  to  get 
across  the  shear  layer  vs.  its  axial 
coordinate  at  t = 0,  for  two  roll  rates. 
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Figure  9.  Trajectories  for  bomblets  la,  lb,  lc,  and  28a  and 
their  velocity  at  t = 120  msec.  The  three  dots 
on  each  trajectory  are  at  t = 80,  100,  and  120 
msec . 
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POSITION  OF  BOMBLET 


Figure  10.  The  solid  curves  through  the  dots  give  the  distribution 
at  t = 100  msec  of  the  bomblets  initially  in  rows  a,  b, 
and  c.  The  significance  of  the  dots  is  given  in  the 
text.  The  dashed  line  is  obtained  from  a test  result. 
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Figure  11.  The  distribution  at  t = 100  msec  of  the  bomblets  initially 
in  rows  a,  b,  and  c for  two  roll  rates.  The  significance 
of  the  dots  is  the  same  as  in  Fig.  9. 
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LIST  OF  SYMBOLS 


a 

initial  radial  coordinate  of  bomblet  center  of  mass  (m) 

d 

m 

diameter  of  bomblet  (m) 

mass  of  bomblet  (kg) 

P 

pressure  at  the  vortex  sheet  (MPa) 

Poo 

free  stream  pressure  (MPa) 

q 

dynamic  pressure  (MPa)  = % p^V^2 

ro 

r-coordinate  of  the  center  of  mass  of  the  bomblet  when 
it  is  tangent  to  the  shear  layer  (m) 

s ,n 

coordinates  perpendicular  and  parallel  to  the  shear 
layer  (m) 

t 

time  (sec,  msec) 

t 

e 

ejection  time  (sec,  msec) 

t . 
i 

time  it  takes  the  bomblet  to  move  up  to  the  shear  layer 
(sec,  msec) 

x,  r,  <f> 

cylindrical  polar  coordinates  fixed  with  respect  to  the 
missile  (m) 

xo 

initial  axial  coordinate  of  the  bomblet  center  of  mass 
(m) 

A 

cross-sectional  area  of  bomblet  or  missile  (m2) 

CD 

drag  coefficient 

C1 

constant  in  equation  preceding  equation  (5)  (=  a2w) 

V Fn 

forces  perpendicular  and  parallel  to  shear  layer  (N) 

F » F u 

sm  nh 

F , respectively  at  that  position  where  the  bomblet 

center  is  on  the  vortex  sheet  (N) 

F , F , F, 
x’  r’  <p 

components  of  applied  force  in  cylindrical  coordinates 
(N) 
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LIST  OF  SYMBOLS  (continued) 


Mach  number  outside  vortex  sheet 

free  stream  (or  missile)  Mach  number 

velocity  external  to  the  vortex  sheet  (m/sec) 

free  stream  velocity  or  missile  velocity  at  event 
(m/sec) 

weight  of  bomblet  or  missile  (N) 

atmospheric  density  at  event  height  (kg/m3) 
roll  rate  (rad/s) 

shear  layer  angle  with  respect  to  the  x-axis 
roll  rate  (rev/sec) 

asymptotic  value  of  roll  rate  (rev/sec) 
differentiation  with  respect  to  time  (d/dt) 
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